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Abstract 

The radial part of Klein-Gordon equation is solved for the Woods-Saxon potential within the 
framework of an approximation to the centrifugal barrier. The bound states and the corresponding 
normalized eigenfunctions of the Woods-Saxon potential are computed by using the Nikiforov- 
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I. INTRODUCTION 



In recent years, the exact or approximate solutions of wave equations have received a 
great attention. The Woods-Saxon (WS) potential given by [1] 

V ^= ~ i + J°-n,)/a » 
is one of the much studied potential in quantum mechanical problems [2]. V is the po- 
tential depth, the parameter a is the thickness of surface, and we denote the width of the 
potential by tq , which is proportional with target mass number A. The coupled-channels 
approach is a powerful tool to reproduce in heavy-ion physics, and Woods-Saxon potential, 
as a internuclear potential, has an important role in the coupled-channels calculations [2]. 
Further, the nuclear optical-model potential is widely used to study elastic scattering pro- 
cesses of nucleons, and heavy particles, and generate distorted waves in nuclear reactions. 
The Woods-Saxon potential is one of the three parts of nucleon-nucleus potential in the view 
of optical-model [3]. 

The wave equations with the Woods-Saxon potential can be solved analytically for s-waves 
due to the centrifugal potential barrier, and these solutions including the wave functions have 
been obtained by using different methods [4, 9]. In this work, we give the energy eigenvalues 
and the corresponding eigenfunctions of the radial Klein-Gordon (KG) equation for usual 
Woods-Saxon potential for any £ values by using an approximate term instead of the cen- 
trifugal potential barrier. We solve the radial part of KG equation by using Nikiforov-Uvarov 
(NU) method, which is a powerful method to solve the second-order, linear differential equa- 
tions [10]. Further, we also study the energy spectrum, and the corresponding wave functions 
of the Schrodinger equation for any £ value by using the same approximation. So, we check 
out the consistency of our new approximation scheme in the non-relativistic case. 

The work is organized as follows. In Section II, we solve the radial part of KG equation for 
usual Woods-Saxon potential for any £ state by using an approximate potential term replaced 
by centrifugal potential barrier. We find out the eigenvalues and corresponding normalized 
eigenfunctions, and also give the results for s-waves. We also apply our approximation to 
the case of the Schrodinger equation, and it makes possible to control the accuracy of our 
approximation in non-relativistic region. We summarize our concluding remarks in Section 
III. 



2 



II. BOUND STATES AND NIKIFOROV-UVAROV METHOD 

In spherical coordinates, the radial part of KG equation can be written as [11] 

r h 2 d 2 h 2 £(£+l) 1 . „ , /ri Tr . ,,„.-, ,. . 

where .E is the energy of the particle, mo is the rest mass of particle, and t is the angular 
momentum quantum number. 
Let us write the potential as 



V(x) = -V , (3) 

V > 1 + e f3x ' V ) 

where x = r — r , and (3 is a short notation, i.e., /3 = l/o. Eq. (2) can not be solved 
exactly because of the centrifugal potential, but the term can be expand about x = as the 
following 

h 2 £(£+l) D x ,x, 9 

2m r 2 (! + ^) r o V y 

because the nuclear distance r can not fluctuate very far from the equilibrium for rather 

high vibrational levels [12, 13, 14], which means that the series expansion in Eq. (4) is 

valid for small x values. The parameter D in the above equation is a short notation, i.e., 

jj _ h 2 £(£+l) 

We prefer the following form instead of the centrifugal potential barrier 



V[{x) 



Do + 



D 1 



+ 



Do 



D 



(5) 



1 + e? x (1 + e^) 2 

where we use three new parameters D , D 1 , and D 2 . The parameter D corresponds to 
the constant term in series expansion in Eq. (4), and the remaining two terms Di,D 2 
correspond to term proportional with 1/(1 + e^ x ) , and proportional with 1/(1 + e@ x ) 2 in Eq. 
(2), respectively. 

Expanding the potential V{(x) around x = under the same condition, and than com- 
bining equal powers with Eq. (4), one can find the arbitrary constants D^i = 0,1, 2) in the 
new form of the potential as following 



12 4 

48 _8_ 

1 ~ /^V 2 + (3r ' 

So, we get two different 'effective' potentials of the form 

In Figs. 1, and 2, we plot the variation of the 'effective' potentials V e ff(x) , and V^f(x) 
with respect to f3x for three different i values. The Figs. 1, and 2 show that there is a well 
consistency between of the potential V e ff(x) , and our approximation in Eq. (5), where we 
use the numerical values V = 43.1 MeV, a p = 0.67 fm, a a = 0.55 fm, r = 3.44731 fm, 
m p = 1.007825 amu, m a = 1.00866 amu for particle, and antiparticle, respectively [19]. We 
see that the deviation from the effective potential V e ff(x) appears for higher £ values (£ > 5) 
for the range starting from /3(r — tq) >~ 2.5. So, if we use the geometric average value 
A = 56 [20], we get an average as r = 4.91623 , and taking a = 0.654 fm [20], we could set 
an upper limit such as r < ~ 1.70974 fm for a good consistency of our approximation. 

Substituting Eq. (5) into Eq. (2), we get 



4- |1M>BV„ - ^DD,]^ + \SX - IpDDAj^l^M = , (9) 
By using the transformation z = 2(1 + e 13 ^ 1 (0<x<oo^0<2;<l),we have 



d 2 <P(z) + 2(1 -z) d<j>(z) + 



dz 2 z{2-z) dz [z{2-z)Y 
where 



2 2 2 2 

— a x z — a 2 z — a 3 



(f>(z) = 0. (10) 
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- a\ = j 2 (5 2 V 2 - 2m DD 2 /h 2 ) , 

-a\ = 4 (^ 2 EV - 2m DD 1 /h 2 ) , 

P 

-a\= -^(5 2 (E 2 -m 2 c 4 )-2m DD /h 2 ). (11) 
To apply the NU-method, we rewrite Eq. (10) in the following form 

^)+|||^) + ^^) = 0, (12) 

where a(z) and a(z) are polynomials with second-degree, at most, and f(z) is a polynomial 
with first-degree. By using the following transformation for the total wave function 



<P(z) = £{z)1>(z) (13) 
we get a hypergeometric type equation 

a(z)^j"(z) + t(zW(z) + \il>(z) = 0, (14) 
where £(z) satisfies the equation 

e(z)/£(z)=*(z)/a(z). (15) 

and the other part, ip(z), is the hypergeometric type function whose polynomial solutions 
are given by 



^)=-^ ) £ i [o n (z)p{z)], (16) 
where b n is a normalization constant, and the weight function p(z) must satisfy the condition 



j- z [a{z)p{z)]=T{z)p{z). (17) 



The function n(z) and the parameter A required for this method are defined as follows 



. , a'(z)—f(z) Lcr'(z) — f(z) sn s , , s 

= 2 V ( 2 } ~ ^ + W ' ( } 

\ = k + it\z) (19) 

The constant is determined by imposing a condition such that the discriminant under the 
square root should be zero. Thus one gets a new eigenvalue equation 

A = \ n = -nr'- ^-^la\{n = 0,1,2,...) (20) 

where 

t(z) = f(z) + 2n(z) . (21) 

and the derivative of t(z) must be negative. 
Comparing Eq. (10) with Eq. (12), we have 

f{z) = 2(1 - z) , a(z) = z{2 - z) , a(z) = -a\z 2 - a\z - a\ (22) 
Substituting this into Eq. (18), we get 

tt(z) = ±y/(al - k)z 2 + {a\ + 2k)z + of. (23) 

The constant k can be determined by the condition that the discriminant of the expression 
under the square root has to be zero 

(a 2 2 + 2kf -Aal(a\-k) = 0. (24) 

The roots of k are fc lj2 = —\cL2~\ a \ Z! r\ a 3 A, where A = \Ja\ + 2a^ + 4af . Substituting 
these values into Eq. (23), we get for n(z) for ki 
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n(z) = T 



V 9 0. J 



z + a 3 



(25) 



and for /c 2 



ir(z) = T 



(4 + f) 



2 2 

Now we find the polynomial t(z) from ir(z) for /c 2 



z - a 3 



(26) 



r(z) = 2 + 2a 3 -2(| + ^ + l)z. 
so its derivative —2^ + ^ + 1 ) is negative. We have from Eq. (19) 



(27) 



and Eq. (20) gives us 



A = - | («2 + °3 + Aa 3 + A + a s) , 



(28) 



A 



° 3 i 1 A i 2 

+ 1 + n — n . 



(29) 



X n = 2n(- + 2 

Substituting the values of the parameters given by Eq. (11), and setting A = A n , one can 
find the energy eigenvalues for any Estate 



2^2 , /1X2t/2 8m D , 



[(3 2 n 2 + 45 2 V< 



D 1 + D 2 )}V 



2{{3 2 n\ + 45 2 V 2 ) 



± 



(3n\ 



2m 2 c 4 S 2 + 2 -^{2D + A + 
2(/3 2 n 2 + 45 2 V 2 ) 



ggCgH-gj 
/3 2 n 2 + 45 2 \/ 2 



16 



, (30) 



where n is the principal quantum number, and 



rai = -(2ra + 1) + + 4a 2 . (31) 

From this result, we can easily get the energy spectra for s-waves by setting D = D = 
D 1 = D 2 = 



7 



E n/=0 ---± — ^ 



rriQC 4 d z 



1 

16 



(32) 



where 



n[ = -(2n + 1) + yjl + 4a? 



(33) 



where a' x = ai(€ = 0). 

It can be seen that the eigenvalues are real under the condition that 

+ . Further, the parameters D 0) Di,D 2 used 



2m 2 c 4 5 2 +^-(2Do+D l+ D 2 ) 



2(^+452^) 



> 



^(Pl+J>2) l- 



to describe the approximate potential form in Eq. (5) are real in the case. 

We give the variation of the bound states energy of particle + a nucleon with A, and 
antiparticle + a nucleon with A with respect to n for different £ values for in Figs. 3, and 
4, respectively. We choose the target mass number as A = 20 , and use the numerical values 
given above. 

In order to find the eigenfunctions, we first compute the weight function from Eq. (17) 



p(z)=z a >(2-z) A , 



(34) 



and the wave function becomes 



z n+a 3 ( 2 _ Z y+A 



(35) 



z a *(2 - z) A dz n 

where b n is a normalization constant. The polynomial solutions can be written in terms of 
the Jacobi polynomials [15-17] 



iP ne (z) = b n Pi a *> A ^(l-z), A>-1, a 3 >-l. 
On the other hand, the other part of the wave function is obtained from Eq. (15) as 



(36) 



£(z) = z a *l 2 (2 - z) A l 2 



(37) 
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Thus, the total eigenfunctions take 



<f> nt (z) = b' n (2 - z) A ' 2 z^ 2 A ^ (1-z). (38) 



where b' n is the new normalization constant. It is obtained from 



To evaluate the integral, we use the following representation of the Jacobi polynomials [16, 
17] 



pM( z)= T ( n + °+ 1 ) T(n+1) Tjn + a + g + m + l) z-1 

n [) nW(n + a + <; + l) ^ T(m + l)T(n - m + I) T(m + a + 1) 1 2 ' ' 

(40) 

Hence, from Eq. (39), and with the help of Eq. (40), we get 



^[g(n, m) x g(r, s)} \b' n \ 2 £ *™+»+«3+i (2 _ z) A-i dz = x ? (41) 

where g(n, m), and g(r, s) are two arbitrary functions of the parameters A, and a 3 , and given 
by 



, (n>m) = ( _ ir+ i/ 2( _ 2) - m+ i/ 2 J^ +w+1 ) 

n!T(A + a 3 + n + 1) 



" r(n + l) r(A + a 3 + n + m+l) 



f T(m+ l)r(n-m+ 1) T(a 3 + A + 1) 



(42) 



and 



g(r, s) = g(n, m) (n — > r; m — > s) . (43) 

The integral in Eq. (41) can be evaluated by using the following integral representation of 
hypergeometric type function 2 F\ (a, b; c; z) [13] 



2 F 1 (a,b;c;z)= T{b ^_ b) £ (1 - t) M (1 - tz)- dt , (44) 



by setting the variable z — > §, and taking c— 1 + b ,z — 1, one gets 



/ o 1 i -(2- f )-"*=™^^(»,M + ^), (45) 



'(1 + 6) 
From last equation 

- 1 ..m+s+a.+i ,„ ... r(m + s + a 3 + 2)r(l) 



/ ^ m+s+a3+1 (2 - ^) A - X rf2 

JO 



o 2r(m + s + a 3 + 3) 



1 



(46) 



x 2 Fi(l - A, m + s + a 3 + 2; m + s + a 3 + 3; - ) , 



where we set b — 1 = m + s + a 3 + 1 , and a = 1 — A. 

By using the following identities of hypergeometric type functions [18] 



r(i& + i)r(&-a + i) 

2 F 1 (a,6;c;-l) = V2 J, -y , (a - 6 + c = 1 , 6 > 0) (47) 

2 Fx (a, 6; c; ^) = T 2 F 1 (a, c - &; c; -1) , (48) 
the hypergeometric type function of 2 F 1 (a, b; c; z) in Eq. (46) can be evaluated as 



2 iq(i-AM + ^)= ^ rf{ + + A }y ( A = b - l )> (49) 

Finally, we get the normalization constant as 



= 2^ r(m + S + a 3 + 3)r(| +A) 

1 nl ~ /3 v ^Fr(m + s + a 3 + 2)r(l + A)[( ? (n,m) x g(r,s)] ' 1 J 

It is interested to study the approximation scheme of potential given by Eq. (5) in the 
case of Schrodinger equation (SE). It also makes to possible to check out the results obtained 
by using the approximation Eq. (5) in the non-relativistic region. 
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The radial part of SE is given by 



(51) 



2m dr 2 [ 2m r 2 , 

By using the approximate expression of the centrifugal potential barrier, and the trans- 
formation z — (1 + e^)" 1 , we have 



d 2 (j)(z) 2z-l d(j)(z) 



dz 2 



z(z-l) dz [z(z- l)p 



2 2 2 2 

— e — 7 z — « z 



^(z)=0. 



(52) 



where 



e = 



2m 



(DDo - £) , 



(3 2 h 2 



2m 



DD 2 . 



(53) 



Following the same procedure, the energy eigenvalues are written 



E n f — 



ti 2 



2m r 2 



-Da 



h 2 



2m a 2 



2n+ 1 + 



\ 



, , m + 1)« 2 „ 

A H o -^2 



r 2 



l(l+l)q 2 / £^ _|_ ^ \ _ 2m Vba 2 



2n+l + Jl + W^D 2 



.(54) 



We plot the variation of the energy eigenvalues obtained from Eq. (54) with respect to 
n for t = in Fig. 5. We give the dependence of the energy spectrum to n for the values 
£ = 1,2 in Fig. 6 by using the numerical values given above. 

We get easily the energy spectrum for s-waves 



E, 



h 2 



nl=0 — 



2m a 2 



m 2 + ( 



moV a 2 
k 2 {n + 1) 



2 m V a 2 

+ h 2 



(55) 



which is consistent with obtained in Ref. [9]. 
The corresponding eigenfunctions are written 
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Mz) = a n (l- z)~ A z-*P^>- 2A) (1 - 2z) . (56) 

where A = ^e 2 + 7 2 + k 2 , and a n is normalization constant. Let us discuss the behavior of 
the approximate wave function at the origin. Taking into account that z — (1 + e@ x )~ l and 
x = r — r , we get the following limits 



if x — > oo 

(57) 

| if a: ^ 

We should take the limit z — > | to discuss the behaviour of the wave function at origin. We 
obtain the wave function in this limit as 



<Pm (z) ~ (i)- A - e Pi- 2e '- 2A) (1 - 2* -> 0) , (58) 
The Jocabi polynomials has the following form in this limit [17,18] 



Pt 2e '- 2A) (y) = 





(n- 


2/1 \ 


I k 1 







;i + y)*( y -i)»-*, (59) 



where y identifies 1 — 2z . So we write the approximate wave function in the first order of z 



where 



4>nt (z) 



{(_!)"-* + [(-l)"-*" 1 ^ — k) + (-l) n - k k\ z^r ... }f(n, k, e, A) , 



(60) 



1 n 
- V 


f n - 2e ^ 




( n - 2A \ 


On 

Z fc=0 


I k J 




\^ n — k J 



(61) 



It is seen from Eq. (60) that the approximate wave function is finite at the origin. It will 
have the same behavior with the approximate solution. 
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As a final test, we compare our results obtained from Eq. (54) with the ones given in 
Ref. [21], where the problem is solved numerically, and the results are obtained by using a 
MATLAB package called as the MATSLISE. So, we use the following form of the Woods- 
Saxon potential in Eq. (3) 

50 

^) = - 1 + e( ,-r)/o.6 . ( 62 ) 

where we should stress that the variable x in last equation denotes the radial coordinate r 
in the rest of the present work. 

We obtain energy eigenvalues of the Woods-Saxon potential given in Eq. (62) by using 
the approach in Eq. (5) numerically, and summarize the results in Table I. We give the 
eigenfunctions obtained without, and with using the approximation given in Eq. (5) for 
the Woods-Saxon potential given in Eq. (62) in Fig. 7 for n — 1 , and I — 1 , and Fig. 8 
for n = 6 , and £ = 5 , respectively. The deviations between the results obtained with and 
without the approximation increase as the values of n , and £ increase, and there is a shifting 
to the left between the approximate and exact eigenfunctions. These are expected results, 
because the validity of the approximation is specified in the Figs. 1, and 2 for the range as 
(r <~ 1.70794 fm), while we set the molecular distance r as r = 7 in Eq. (62). 

III. CONCLUSION 

We have solved analytically the radial part of the KG equation for the usual Woods- 
Saxon potential in the framework of an approximation to the centrifugal potential term 
for any £ values. The energy spectra and the corresponding wave functions are obtained 
by applying the NU- method. We give in Figs. 1, and 2 the variations of V e ff(r) , and 
K//( r ) w ^ n respect to f3(r — r ) for £ = 1,2,5. We have pointed out that the consistency 
between V e ff , and our new approximation is very well, and reliable results can be obtained 
by using the approximation scheme used in the present work. To check our results, we have 
also calculated analytically the energy eigenvalues of the particle and antiparticles for the 
s-waves. We have found that the analytical results are consistent with those in Ref. [4] if 
q — > 1. Further, we have also studied independently the energy spectrum of the Schrodinger 
equation by using the same approximation to check out our results in the non-relativistic 
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region. We have seen that the results obtained analytically for the case of the Schrodinger 
equation are the same for s-waves with Ref. [9]. Further, we give some numerical results for 
the eigenvalues of the Schrodinger equation with different values of the quantum numbers 
n , and t in Table I. We compare the eigenfunctions for two different I values obtained with, 
and without using the new approximation in Figs. 7-8. 
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FIG. 3: The dependence of a particle energy levels in the case of the Klein-Gordon to n for different 
I values 1,2,3,4. 
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FIG. 4: The dependence of a antiparticle energy levels in the case of the Klein-Gordon to n for 
different i values 1, 2, 3, 4. 
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FIG. 5: The variation of the Schrodinger energy level with respect to n for i = . 




n 



FIG. 6: The variation of the Schrodinger energy levels with respect to n for t = 1, 2. 
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FIG. 7: The Schrodinger eigenfunctions without (full line), and with using the approximation in 
Eq. (5) for Woods-Saxon potential given in Eq. (57) for n = 1 , and t = 1 . 




